Non-Gaussianities in models with a varying inflaton decay rate 
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We consider the expected level of primordial non-Gaussianities in models in which density per- 
turbations are produced by spatial fluctuations in the decay rate of the inflaton. We consider 
both the non-Gaussianities resulting from the self-couplings of the field that controls the decay 
rate as well as from the non-linear relation between field and curvature perturbations. We show 
that in these scenario non-Gaussianities are of the "local" form, ie well described be the ansatz 
1Z — lZ g + /jviC^-g — (TZg)). This is a consequence of the fact that they were created when modes 
were already outside the horizon. We show that f§f L is naturally of order a few in these models, 



much larger than what is expected in the standard one field models of inflation (/ 
possibly accessible to observations. 
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I. INTRODUCTION 

It is safe to say that the era of "precision cosmology" is 
upon us. The latest CMB anisotropy measurements from 
WMAP [1] have underscored how successful our standard 
cosmological model is. Many of the parameters describ- 
ing the matter content of the universe are now firmly es- 
tablished [2] . Moreover the conclusion that the structure 
we observe in the universe grew from seeds that were in 
place much before recombination seems almost unavoid- 
able [3]. 

The leading scenario for explaining how these initial 
seeds were created is inflation. In this picture the uni- 
verse went through a period of accelerated expansion very 
early on. During inflation the energy density is domi- 
nated by the potential energy of a slowly rolling scalar 
field, the inflaton. During this period quantum fluctua- 
tions in the inflaton field were stretched outside the hori- 
zon and became the seeds for structure. The predictions 
from this scenario have been worked out to great detail 
in the last few decades (see for example [4,5] and refer- 
ences therein). Basically most models of inflation pre- 
dict a scale invariant spectrum of Gaussian fluctuations 
of the curvature. A stochastic background of gravita- 
tional waves is also predicted, although its level is quite 
uncertain. 

The Gaussianity of the perturbations in the standard 
picture are a consequence of the slow-roll conditions on 
the inflaton potential. These conditions are required for 
the field to be slowly rolling and the universe to be inflat- 
ing. The required flatness of the potential implies that 
non-Gaussianities are very small, probably below what 
could be observable (see [6] and references therein). If 
the inflaton field has several components it is possible to 
construct models with larger non-Gaussianities even in 
the adiabatic component. Basically large couplings are 
allowed for the "isocurvature" directions which if suit- 
ably coupled to the adiabatic direction could lead to non- 
Gaussianities in that direction as well (eg. [7-9]). 

There are basically two modifications of the standard 
inflationary scenario in which the curvature perturba- 
tions we observe today were not created during inflaton 



but later. These are the "curvaton" scenario [10-15] and 
the scenario proposed in [16] (see also [17,18]). Both re- 
quire the existence of other light fields during inflation. 
In the "curvaton" model one of these fields eventually 
comes to dominate or at least be a significant contribu- 
tor to the energy density and reheats the universe again. 
The curvature perturbations we observe are a result of 
the fluctuations in the curvaton field as opposed to the 
inflaton. Fluctuations in the curvaton were still produced 
during inflation by the standard mechanism but they get 
converted into curvature perturbations when the curva- 
ton becomes dominant. 

The scenario presented in [16] is different. At the end 
of inflation the energy stored in the inflaton has to be 
converted into normal particles, reheating the universe 
and starting the standard phase of the hot big bang. In 
[16] we suggested that if the efficiency with which the 
inflaton reheats the universe, its decay rate T varied in 
space, density perturbations would be generated during 
reheating independently of those generated by the stan- 
dard inflationary mechanism. If two different regions of 
the universe had different Ts then effectively the inflaton 
would decay into radiation first in one region and then in 
the other. During the time one region is filled with radia- 
tion while the other one is not, the universe expands at a 
different rate in each region, resulting in density pertur- 
bations when reheating is finished. In a similar spirit, the 
in-homogeneous freeze out of an unstable particle whose 
mass and decay rate are determined by a light scalar field 
and that briefly dominates the energy density would also 
lead to density perturbations [17]. 

The motivation for these new scenarios is that in gen- 
eral the decay rate of the inflaton as well as the masses 
and decay rates of particles are determined by the expec- 
tation values of scalar fields. If those scalar fields were 
light during inflation they fluctuated, which lead to den- 
sity perturbations. 

In this paper we will analyze under what conditions 
the fluctuations in the scenario of [16] are Gaussian. For 
example, the potential for the fields that determine the 
inflaton decay rate do not have to satisfy any slow-roll 
conditions so one might expect their fluctuations to be 
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significantly non-Gaussian. We will find that gcncrically 
we predict non-Gaussianities that are larger than those 
predicted by the standard inflationary model although 
within the observational bounds. Moreover these addi- 
tional non-Gaussianities could be accesible to future ex- 
periments. There are two avenues for searching for pri- 
mordial non-Gaussianities observationally, studies of the 
CMB [19] and studies of Large Scale structure as traced 
say by galaxies or weak lensing (sec [20] for a comprehen- 
sive review of non-Gaussinanity in the context of Large 
Scale Structure). 

The paper is organized as follows. In §11 we summarize 
how the curvature perturbations are generated in this 
scenario. In §111 we discuss the different sources of non- 
Gaussianities and finally in §IV we make contact with 
observations. 



II. DENSITY PERTURBATIONS IN THE NEW 
SCENARIO 

The decay rate of the inflaton (which we call S) can 
be written schematically as, 



r = m s g 2 K, 



(1) 



where ms is the mass of the inflaton, g its coupling con- 
stant and K are kinematic suppression factors of the form 
(1 — m d /m 2 s ) where nid stands for the masses of the de- 
cay products which are also functions of the scalar fields. 
All three quantities, ms, g and K are functions of the 
scalar fields in the theory so clearly there are a wide rage 
of possible specific models in which V has a different de- 
pendance on the fields 

During inflation the light fields </>j have expectation 
values 4>i around which they fluctuate. To keep our dis- 
cussion generic we will express the decay rate as an ex- 
pansion in powers of those fluctuations, 



dv !± i 9 2 r 



biS(f>j + ■ 



(2) 



In [16] it was shown that to linear order in the de- 
cay rate fluctuations ST/T curvature fluctuations 1Z are 
generated. We found, 



(3) 



where the constant a depends on the ratio of f to the 
Hubble expansion at the end of inflation (H). As T/H 
increases a decreases. It is clear for example that in 
the limit T/H » 1 reheating happens almost instanta- 
neously in all regions regardless of the value of T and thus 
the resulting density perturbations will be very small. In 
[16] it was proved that in the limit when L << H 



(4) 



We will now present a different derivation of the same 
result that will allow us to extend our treatment to in- 
vestigate non-Gaussianities in a later section. 

We first consider reheating in a homogeneous universe, 
ie with a constant T. We will treat the inflation as non- 
rclativistic matter when it oscillates around the minimum 
of its potential after inflation ends. We will call e m the 
energy density in the oscillating inflaton. During reheat- 
ing, the energy density in the inflaton is converted into 
radiation (with energy density e r ). The background met- 
ric is that of a flat FRW universe, ds 2 = -dt 2 + a 2 (t)dx 2 , 
a is the scale factor. The Hubble constant is H = a/ a, 
with ' = d/dt where t is the time. The evolution equations 
are 



-ZHe r . 



H 2 = 



-AHe r + Te m , 
8ttG 



(5) 



After inflation the fraction x = e m /(e m + e r ) = 1, and 
remains constant until it suddenly drops to zero when 
H - T. After that point a oc t 1 / 2 . 

The key point to understand is that universes with 
different values of T will have expanded by different 
amounts by the time reheating in over. In figure 1 we 
show the final scale factor as a function of F. We plot 
ln(#) = ^(at- 1 / 2 ) vs. \n(Y/H) where H is the expansion 
rate at the end of inflation. The asymptote at T/H — > 
corresponds to g oc T -1 / 6 , a result that can be obtained 
analytically [21]. 
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FIG. 1. Expansion factor normalization g as a function of Y / H 
after reheating. The curve was obtained by solving equations (5) 
numerically starting with e r = and Hubble expansion H at the 
initial time. 
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After reheating is over the metric is ds 2 — —dt 2 + 
g 2 (T)tdx 2 where <? 2 (r) is the function plotted in figure 1. 
To obtain the density perturbations in our model we just 
need to promote T to a function of position, 



n = n gl + f^i 1 {n 2 -{n 2 )) 



ds 2 = -dt 2 + g 2 {T{x))tdx 2 



+ fN 9 L (K2 <^2» + /£r 92 ft»lft»2 + • • • (9) 

where 7Zi g and 7Z2 g are Gaussian random fields and 
(6) the fj$ L control the level of non-Gaussianities. We 



This metric will be a solution of Einstein's equations if V 
is sufficiently slowly varying, that is to say if T changes 
on scales much larger than the horizon. Clearly this is a 
very good approximation as we are interested in changes 
of T over scales that are of cosmological interest today. 
It is important to realize that we have not made an ex- 
pansion in powers of ST/T, we have just assumed that 
the gradients are small. As a result we can use the above 
metric to calculate quantities to higher order in 8T/T as 
long as we do so when the variations of T are on super- 
horizon scales. 

To make contact with observations we need to calculate 
the spatial curvature of surfaces of constant time. For the 
metric in equation (6) we obtain [5,23], 

"«TT- ( ? ) 



1Z= — 



where a = — dlng/dlna is just minus the slope of the 
curve in figure 1 at the appropriate value of T. This 
result was already obtained in [16] by the standard per- 
turbation theory methods. In linear perturbation theory 
the variable 1Z is related to the curvature of hypersurfaces 
of constant t (R^) by (see for example [5]) 



i?(3) = -4v 2 ^. 



(8) 



We will generalize this expression to include higher order 
terms later. 

A few things are worth noting. If after reheating the 
composition is the same everywhere, if we are dealing 
with adiabatic perturbations, the universe will expand 
by the same amount everywhere. As a result ratios of 
the scale factor between different places are constant in 
time, and thus 1Z is constant. Certainly because the cou- 
plings and masses of particles are different in different 
places at reheating in this scenario there is the possibility 
of generating different ratios of abundances of particles 
as a function of location. Thus our scenario can cre- 
ate "isocurvature" density perturbations as well. This is 
similar to what happens in the "curvaton" scenario. Al- 
though this is a possibility, it is rather unlikely. Isocurva- 
ture perturbations could result if either there is varying 
decay into species out of equilibrium or if the decay of the 
inflaton lead to baryon asymmetry or lepton asymmetry 
fluctuations. 



have departed slightly from the notation of Komatsu & 
Spergel [19] by working with the curvature 1Z rather than 
the gravitational potential because 1Z is constant out- 
side the horizon while the gravitational potential evolves 
whenever the effective equation of state of the universe 
changes. We have also allowed for the possibility that 
more than one Gaussian field TZ\ g and lZ2 g enter into the 
quadratic corrections. As equation (2) makes clear, this 
could happen if more than one field <j>i determined the 
decay rate of the inflaton. The first thing to point out is 
that the extra terms do not enter into the calculation of 
the three point function to lowest order in $nl so obscr- 
vationally they are not relevant. Higher order terms are 
suppressed by extra powers of 1Z <~ 10~ 5 . 

There are several sources for the non-Gaussianities in 
1Z. The fluctuations 8cj>i in equation (2) could themselves 
be non-Gaussian. Even if those were Gaussian the rela- 
tion between 8<pi and the decay rate T in (2) is non-linear. 
The fact that T is by definition positive means that the 
relation cannot be purely linear, the size of those non- 
linearities however varies from model to model. More- 
over as shown in the previous section g(T) is a non-linear 
function of T which leads to additional non-Gaussianities. 
We will consider the different cases in turn. 



A. Gaussian 8(f>i fluctuations 

We start by considering a simple example in which 
there is a field <f> that sets the coupling of one of the 
channels through which the inflation can decay and that 
in addition there is an extra channel that does not fluc- 
tuate. This example was considered in [16]. The decay 
rate can be written as, 



so that 



r = r + r 1 (5) 2 

= T 0+ T l{ l + 5 4) 2 , 



r-f = ri u 8<t> 2 - (8<t> 2 ) 
r r 1 </> 4> 2 ' 

f = r + r l( i + 



(10) 



(ii) 



III. NON-GAUSSIANITIES 

In this section we will discuss the level of non- 
Gaussianities that we can expect in this scenario. We 
will consider the curvature perturbations 1Z and write, 



The ratio Ti/T determines the fraction of the decay rate 
that is controlled by <f). 

We first consider the transformation between T and 1Z 
only to linear order, 



3 



r v 



+ 



Sep 2 - (S(j) 2 } , 



(12) 



where a is less than 1/6 (see figure 1). We can rewrite 
equation (12) as 



n = iz g + f% L (iz 2 q -(iz 2 q )) 



f. 



n 
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(13) 



Equation (13) shows that the smaller the efficiency of 
our mechanism of generating density perturbations, the 
larger the f^ L - The inefficiency in the present example 
comes from two sources, a small value of a which occurs 
when the decay rate of the inflaton in not much smaller 
than the Hubble constant at the end of inflation. The 
other source of inefficiency comes form the fact that the 
field only controls part of the channels through which 
the inflaton can decay. Again, the smaller the fraction 
of the decay rate controlled by 0, the largest the non- 
Gaussianities. 

The origin of the relation between inefficiency and non- 
Gaussinanitics is simple to understand. If we keep the 
level of fluctuations in 1Z fixed, which me must to ex- 
plain the observed level of density perturbations, as we 
decrease the efficiency we must increase the level of <50/0. 
As a result quadratic terms in 8<f>/(f> become more im- 
portant relative to the linear term and thus the non- 
Gaussianity increases. As this arguments makes clear, 
the relation between inefficiency and non-Gaussianity is 
more general than the particular example we have chosen 
for the form of the decay rate. 

We now need to consider the fact that the relation be- 
tween the decay rate and the final curvature is non-linear, 
that is the final value of the scale factor is a non-linear 
function of the T as figure 1 shows. To be more spe- 
cific we will calculate the curvature of the hypersurfaces 
of constant time in the metric of equation (6). We can 
do so without doing a linear approximation in (5F/F al- 
though we are calculating only the leading contribution 
in an expansion in terms of the horizon scale over the 
scale of variation of T. We get, 



r [V 2 ln 5 -||Vln 5 || 2 ], 



(14) 



where a = g(x)t 1 ^ 2 . The first term is the direct gener- 
alization of the linear perturbation formula, with 1Z = 
\n{g/g). As it should be only spatial changes the loga- 
rithm of g signal curvature perturbation. 

We are interested in calculating the curvature to sec- 
ond order in ST/T. There will be two contributions, one 
from expanding the Ing to second order and the other 
from the second term in equation (14). These second 
term arises because General Relativity itself is a non- 
linear theory. We will comment on these shortly. For the 
moment we focus on the first term. With the definition 
1Z = \n(g/g) we have, 



K=-a 



ST a + f3 5T 2 - (ST 2 ) 



+ ■ 



(15) 



where a = — (ilna/dlnT, the slope of the curve in figure 
1 and (3 is the curvature (3 = d 2 In a/d In Y 2 . Thus the 
fact that the relation between T and total expansion is 
non-linear adds an additional contribution to f^ L , 
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a + P 
2a 2 



(16) 



This contribution is at least equal to f^ L — 3 and in- 
creases as the process becomes more inefficient (a < 1/6) 
because <5r/F increases to produce the same level of 
curvature fluctuations. Note however that these non- 
Gaussianities do not include the extra enhancement fac- 
tor T/T±. The point is that the inefficiency in the transla- 
tion between ST/T and 50/0 makes 50/0 larger and thus 
terms quadratic in 50/0 lead to larger non-Gaussianitics 
than terms arising from 5r/F corrections. 

We now comment on the second term in equation 
(14). It contributes an additional effect that leads to 
an fjf L r~j — 1, thus it is subdominant to what comes 
from the first term and is independent of a. Moreover as 
the modes enter the horizon additional non-Gaussianities 
from non-linearities caused by term we have neglected in 
Einstein's equations will develop. Perhaps more impor- 
tantly one should solve for the observables (such as the 
CMB anisotropies) to second order in 1Z to catch all the 
terms that contribute at the level of f^ L <~ 1. R^ may 
not be the quantity directly relevant to observations. It is 
beyond the scope of this paper to characterize all these 
terms or to extend the usual linear treatment of CMB 
anisotropies to include them. What we have calculated 
should be interpreted as corrections in addition to what 
GR already gives. 

Finally we want to comment on the possibility that the 
quadratic corrections to the Gaussian curvature are not 
proportional to the linear contribution squared. In our 
scenario this happens trivially when more than one light 
field determines T. For example if 



r = r + r 1 (fi) 2 + r 2 (f 2 -) 2 



(17) 



then there are two different Gaussian fields that enter in 
the quadratic corrections, 

2Ti S(f>i 2T 2 502 

1 01 1 02 



1Z 



<;2 



2r 2 (50i 2ri . 



(18) 



The amplitudes of the non-linear terms defined in equa- 
tion (9) are, 

» [(r,/f)3 + ( r 2 /f fl^ (sin 3 9 + cos39) 
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, Kg2 [(iyiy + tiyf) 2 ]- 1 / 2 

= - sin 6* cos 6>(sm 6* + cos 0) 
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[(Tx/f^ + ^/f) 



21-1/2 



■ sin#cos#(sin# — cos 6*), 



with cos# = iytrf + r?] 1 /* an d sine = iyprj + r 2 ] 1 / 2 . 

Other than the slight redefinition of /^£ x these addi- 
tional terms are unimportant if we are only interested 
in calculating effects to first order in the /atz,s. The 
other terms will only enter at higher orders and will be 
thus suppressed by extra powers of 1Z ~ 10 . Gener- 
ically then, modifications of the standard Jnl form, 
1Z = lZ g + /nlC^ ~ are what is expected in our 

scenario. It turns out this remains true even when self 
couplings of the field are considered. 



B. non-Gaussian 



fluctuations 



g N (l,---N) = -i J dt(0\{H(l)5ct>(2)---Sct>(N),H mt }\0) 

The N-point function is most easily expressed in 
Fourier space. For a massless field in de-Sitter we ob- 
tain, 

g(k u ■ ■ •k JV )= (2^) 3 «5 D (k 1 + • • • + kjv) x 
H 2N k 3 

VW^£rI N {k 1 ,k 2 ,---k N ) (24) 



I v = / 2Re[-*(l - i*U7) ■■■(!- iftjvijje'*"'] 



with fct = fci + fc2 + - • - + kN and r/ re /i is the conformal time 
of rchating. Following [6] to project onto the vacuum of 
the interacting theory we deform the time integration so 
it has some evolution in cuclidcan time. We explicitly 
give the results for the 3- and 4-point functions, 



In this section we calculate the non-Gaussianities in 
the S(f> fluctuations if the potential for <f> is not quadratic. 
We will follow [6] , although our situation is much simpler 
because the field makes a negligible contribution to the 
energy density so we do not need to take into account 
perturbations in the metric . The action for 5<j) is, 



= J a 4 d 3 xdi] 

V m y(3) 



- (V<50) 2 



2a 2 



2! 



3! 



IT 



(20) 



We denote the fcth derivative of the potential. 
We quantize the field in the usual way to obtain, 



d 3 k 

(2tt)3 



[/(fc,77) e -*- x a k + r(M)e ik ' x 4] (21) 



where [aku^ ] = (2ir) 3 S D (k x — k 2 ) and f(k,r]) are the 
mode functions. For the massless field in de-Sitter space 
(V^ = and H constant), 



H 



:(1 — ikrj)e 



ikrj 



We 



to calculate the N-point function gjv(l,2, 
points at the time of rchating, 



(22) 

want 
■ N) for N 



j 8 E»<j Mjkj Uiki 

I 4 = 7: — r 



kt 



3^ + Nk ^ 



(25) 



where 7 = 0.577216 is the Eulcr's constant and N kt — 
ln(\k t r] re h\) is the number of e-foldings from the time k t 
crossed the horizon (fc t r; = —1) to the time of reheating. 

The value of In is dominated by the term proportional 
to Nk t because the modes that we observe left the hori- 
zon roughly 60 e-foldings before the end of inflation. This 
has a very simple interpretation. All terms in equation 
(25) are of order one except for that proportional from 
Nk f ■ They represent the effect of the self-coupling on the 
modes as they are being generated, at horizon crossing 
during inflation. The terms proportional to N kt however 
originates from the classical evolution of the field outside 
the horizon. To see this in a concrete example lets con- 
sider the corrections induced by classical outside the hori- 
zon evolution of a field with a small third and fourth order 

coupling. We will write S(f)(x) = <J0( O ) (x) + 5(j>m (x) H 

where <$0( O )( X ) denotes the Gaussian fluctuations that 
are produced in the absence of these couplings, 5(f)(i)(x) 
the first order corrections, etc. We can neglect spatial 
derivatives and assume the field is slowly rolling so that, 



^(3) 

~2T 



~3T 



8f 



(0)- 



(26) 



g N (l, 2, • • • N) = (01U- 1 60(1)6^(2) ■ ■ ■ 5cj>(N)U\0) 



U 



/•treh 

= cxp [—i I 

Jt 



Hint dt] 



(23) 



where t re h is the time of rehating, to is some early initial 

3! 



time and Hint = ^-5(p 6 + ^-5(p + • • • is the Hamilto- 
nian in the interaction picture. To lowest order in H int 
this gives, 



Thus in de-Sitter space (H constant) 



N r^ (3) ^ 2 y {i K^ 1 

(o)+— ^(o)]- 



3iJ 2L 2! 



(27) 



where we have introduced the number of e-foldings N = 
HAt. The three and four point functions implied by 
these expressions are exactly the same as those corre- 
sponding to the terms proportional to N kt in equation 
(25). A discussion of the non-Gaussianites generated by 
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outside the horizon classical evolution of a scalar field, in 
the context of a multicomponcnt inflationary model can 
be found in [8]. 

The physics is very simple, the self-couplings of the 
field creates higher order correlations in two different 
ways. It modifies the the statistics of the fluctuations 
themselves but it also changes the Gaussian part of the 
fluctuations as the field evolves outside the horizon. Be- 
cause we are considering the evolution outside the hori- 
zon it is local in real space, meaning that it leads to 
non-Gaussianitics of the form, 



5(f) = (50(o) + 02^(0) + 03^(0) + 



(28) 



There are deviations from this "local" approximation 
coming from the fact that N kt varies across the scales 
we observed, it changes by roughly 10 % between the 
CMB scales and smaller scales that one could access by 
studying large scale structure. There are also corrections 
of order 1/N kt due to the non-Gaussianities produced at 
the time the fluctuations crossed the horizon during in- 
flation. Those have different configurational depcndance. 

We should also point out that the non-linearities keep 
acting all the way until reheating, even after the universe 
stopped inflating. During the period of matter domina- 
tion, what determines the growth is, 



dt 
H 



2m 



3 

2T 2 ' 



(29) 



where we have just taken H ~ 2/3t during the pe- 
riod between the end of inflation and reheating and 
used H re h = r. Thus we should replace the number 
of e-foldings in equation (27) by an effective number 
N kt = N kt + 3/2(H/T) 2 with H the Hubble constant 
during inflation. 

We will concentrate on the three point function be- 
cause unless (f> is fluctuating very close to the origin, so 
that 5(f)/ <j> ~ 1 the four point function will be smaller. 
In fact one naturally expects (ft 4 ) / '(ft 3 ) (ft 2 ) 1 / 2 - 5<p/4>. 
For our simple one field example the contribution to f^ L 
due to non-linearities in the potential is, 



f n - - 

■I NL — 



r v {3) K~r 

T ia 3H 2 kt 



(30) 



To make an estimate of the size of these non-Gaussian 
contribution we will relate to = m 2 and rewrite, 



Jnl 



J_ V^l 2_rr^ K 
T ia ( m 2 > [ 3H 2 ' 2 • {61 > 



We now notice the in the presence of a mass term the 
evolution of the modes outside the horizon (3HS<fr^ — 
—m 2 5(j) 2 j) creates as scale dependance of the power spec- 
trum because larger wavelengths that came out the hori- 
zon earlier had a longer time to evolve and thus have a 
smaller amplitude at reheating. This is exactly the same 



as what happens in the curvaton scenario [12]. The power 
spectrum spectrum gets a blue tilt from this effect, 



so that, 



rK 

Jnl 



2 m 2 

[ ' ~ 3 H 2 



r v^4> 

— (n-iy 

Via m l 



2 



(32) 



(33) 



Departures from scale invariance are observed to be 
small, that is n — 1 should be smaller than few x 10~ 2 
while Nk t 1S °f order 60-100. If the factor m 2 /3H 2 x 
N kt - 1 and V^^/m 2 - 1 then the non-Gaussianities 
produced by the self couplings are of the same order 
than those produced by the non-linear relation between 
the fluctuating field and the decay rate. To put it an- 
other way, even though V{<j>) does not have to satisfy 
slow- roll conditions, if the spectrum of perturbations in 
nearly scale invariant then the non-Gaussainities are not 
expected to be larger than current observational con- 
straints. 

We can consider a simple example V((f>) — A</> 4 /4! and 
assume that inflation has lasted long enough before the 
scale of our present horizon crossed the horizon during 
inflation, that the field <fr has rolled to the point in the po- 
tential where the change during one e-folding in the clas- 
sical trajectory equals the quantum fluctuations H/2n. 
This point corresponds to ~ \^ X ^H. In that case 
(n — l) mass ~ A 1 / 3 . Thus constraints on the spectral in- 
dex already put A < 10~ 6 . Moreover S(f>/(t> - A 1/3 /27r 
so to match the observed ft <~ 4 x 10 -5 we need 
Tia/f <~ 10~ 4 /A 1//3 . The non-Gaussianity parameter is 
then, 



JN 



NL 



-3 10 3 A 2 / 3 iV, 



kt- 



(34) 



So for example if A 1/3 — 10 2 and N kt <~ 100 we get 



TZ 
NL 



-30 very close to the limit set by WMAP. 



IV. SUMMARY 



We have discussed the different sources of non- 
Gaussianities in the scenario introduced in [16]. In this 
section we summarize the results and make contact with 
observations. In most of the observational studies so far, 
in particular in the latest constraints form WMAP [22] 
the analysis was done in terms of the non-Gaussianities in 
the gravitational potential during the matter dominated 
era, $ = — 3/5ft. Constraints were set on /* L defined 

by, 



* = * 9 + /£l(*2-<*2» + " 



(35) 



Thus f% L = -5/3f% L . WMAP constrained -58 < 
f% L < 134 at 95 % confidence [22]. 
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Non-Gaussainitics in this scenario can have several 
sources. There are contributions that are proportional 
to T/Tict. These come from corrections that are di- 
rectly proportional to 5<p/<f>, either from non-linearities 
in the relation between ST/T and 54>/4> or intrinsic non- 
Gaussianities due to a cubic self interaction of the field. 
In the simple examples we have presented, the first con- 
tribution is f% L = 5/12 (f/IJia). The WMAP con- 
straints already imply that Ti/T > 0.1. The contribution 
from the self couplings, equation (30), combined with the 
WMAP constraints imply that N kt V^/H < 0.1. 

There are non-Gaussianities coming from the non- 
linear relation between TZ and ST/T. These grow as the 
generation of density perturbations becomes more inef- 
ficient, that is as T/H grows but they are not propor- 
tional to T/Y\. In the limit f << H this contribution is 
/* L = —5 and is larger if a < 1/6. 

Finally, all these contributions should be interpreted 
as additions to the f% L <~ 1 coming from second order 
GR effects which we have not considered here. For effects 
produced when modes are crossing the horizon (either on 
the way out or in) spatial gradients will not be negligi- 
ble so the configuration dependance of the three point 
function will be different than in the simple /nl model. 
These corrections are independent of the efficiency of our 
the mechanism to translate 5(f>/(f> into TZ and thus are ex- 
pected to be the smallest. 

A few general conclusions of our study are worth em- 
phasizing. First as noted above, inefficiencies in the 
translation between field fluctuations and curvature per- 
turbations increase the expected non-Gaussianities. Sec- 
ond, the non-Gaussianities in our model are mainly gen- 
erated when the modes of interest are outside the hori- 
zon. As a result it is natural to look for departures from 
Gaussianity that are "local" in real space, meaning of the 
form TZ = 1Z g + ffi L (TZ 2 g - (1l 2 g )) + ■ ■ ■ . If the decay rate is 
determined by more that one fluctuating field there could 
be additional corrections at the second order in TZ, but 
these are irrelevant if one is interested only in N-point 
functions to leading order in /jvl (the actual expansion 
parameter is f^ L TZ so higher order terms are very small) . 
The "local " approximation should be excellent in these 
models. Third because our mechanism is in some sense 
somewhat inefficient (a < 1/6, Ti/T < 1), Jnl ~ few 
is guaranteed in our scenario even if the fluctuations in 
the fields are intrinsically Gaussian. This last point is 
of particular interest because Jnl ~ few corresponds to 
the expected sensitivity of future CMB experimets [19] 
and the standard inflationary scenario predicts correc- 
tions much smaller than these [6] . If the efficiency of the 
mechanism is decreased or if self interactions are impor- 
tant, the level of non-Gaussianities can be increased to 
saturate the current upper limits. However the fact that 
the spectrum of perturbations is close to scale invariant, 
implies that even in the presence of self-couplings the 
natural level of non-Gaussianities is somewhat below the 
current limit although not far from it. 
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